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INTRODUCTION 


Hypothetical Prismatic Channel 

Computational ditficulties associated with the complex valley geometry 
was the dominant consideration in the Teton case study. Irregular valley 
geometry was also an important complicating factor in the Laurel Run case 
study. The hypothetical prismatic channel case study was developed to test 
the model under conditions for which irregular valley geometry was not a major 
concern, The reservoir data for the Teton Reservoir was combined with a 
prismatic channel. The prismatic channel is essentially an extension of the 
S-mile long Teton canyon to 50 miles. 

The prismatic channel, as shown in Figure 1, consists of two reaches of 
constant cross-section joined by a transition reach. The channel cross-~ 
section from mile 0 to mile 5 is constant, and the reach between miles 10 and 
50 has a Slightly wider constant cross-section. Miles 5 to 10 provide the 
transition between the two sections. The cross-section for miles 0 to 5 
approximates the geometry of the narrow, steep-walled Teton canyon just below 
the dam, However, the remaining 45 miles of slightly wider prismatic channel 
does not approximate the wide, flat, abruptly changing topography ot the 
valley of the Teton and Snake Rivers below the Teton canyon. The prismatic 
channel has a constant bottom slope of ten feet per mile and a constant 
Manning roughness coefficient of 0.040 which are representative of the canyon 


of the Teton River. 


Overview of Analysis 
DAMBRK, FLOW SIM 1, FLOW SIM 2, HEC-1, SMPDBK, TRO66, and the HEC 


dimensionless graph procedure were tested using the hypothetical prismatic 


channel. For purposes of comparing models, an initial "base run" set of 
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breach parameters was defined as follows: time to maximum breach size of 1.0 
hour, breach width of 500 feet; and vertical side slopes. However, the HEC 
dimensionless graph procedure is based on an instantaneous compiete removal of 
the dam, and TR66 uses a reservoir water depth versus maximum breach discharge 
relationship based on data from actual past dam failures. The other models 
have breach simulation routines which can reflect the above listed breach 
parameters. However, a solution could not be obtained for these parameter 
values with FLOW SIM 1 and 2 due to instabilities in the computations. 

Computational instability was a major problem in applying FLOW SIM 1 and 
2 and became the prima f--us of the analysis with these models. Obtaining 
solutions was not a p: h DAMBRK. The DAMBRK analysis was used to test 
the sensitivity of the .o wave to variations in breach size and manning 
model results to the parameters NMIN and NSTPS, The TR66 procedure was 
applied alternatively with the peak discharge at the dam computed with the 
TR66 method and DAMBRK, with valley routing by the TR66 procedure in both 
cases. SMPDBK and the dimensionless graph procedure were applied without 


roughness coefficient. The HEC-i analysis included testing the sensitivity of 

particular difficulties being encountered or sensitivity studies being made. 
{ 
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Figure 1 - Channel Cross-Section 
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DAMBRK RESULTS 


The results of the DAMBRK analysis are summarized in Tables 1 throvgn 4 
and Figures 2 through 4, DAMBRK was used to test the sensitivity of the dam 
breach flood wave characteristics to breach width and Manning roughness 
coefficient. A constant manning n is assumed all along the channel. The 
results of four runs are presented in the tables and figures with alternative 
n values of 0.04, 0.06, 0.08, and 0.12 and all other input data constant. 
Likewise, alternative runs are presented with alternative final breach widths 
of 100 feet, 300 feet, 500 feet, and 700 feet and all other input data 
constant. Appendix A is a listing of the DAMBRK base run (breach width of 500 
feet andn values i 0.40). 

Several unsuccessful trial runs were required to debug the input data and 
establish workable distance steps. After these initial runs to get the nodel 
running, computational problems, such as nonconvergence and instability, were 


not encountered, 


doll dl len 


fer 


cl 


Table 1 
DAMBRK - Prismatic 
Peak Discharges 


Maximum Discharge in 1,000 cfs 


perme n. | Lame med “| | aaa als. - i a eeeeei, 4 of ty ors 


Distance : Manning Roughness Coefficient (n) 

From > 0.04 : 0.04 :; 0.04 : 0.04 : 0.06 : 9.08 : 0.12 

Dan : Final Breach Width in Feet 

(miles) : 100 : 300 : 500 : 700 : 500 : 500 : 500 
0 1,224 2,995 3,841 4,563 3,717 3,380 3,108 
5 1,110 2,603 3,468 4,145 3,218 2,838 2,219 
10 1,060 2,349 3,220 3,840 2,840 2,400 1,782 f 
16 995 2,066 2,864 3,374 2,378 1,926 1,368 
22 933 1,897 2,536 2,376 1,978 1,542 1,067 
30 867 1,716 2,135 2,317 1,598 1,241 839 
40 804 1,508 1,777 1,905 1,313 1,027 683 
50 3a 1,357 1,567 1,674 1,136 892 587 


ee See ee ee Se ue 2 8 oe oP SW 


Table 2 
DAMBRK - Prismatic 
Peak Water Surface Elevations 


Maximum Water Surface Elevation in Feet MSL 


Distance : Manning Roughness Coefficient (n) 
From : 9.04 : 0.04 : 0.04 : 0,04 : 0.06 : 0.08 : 0.012 
Dam : Final Breach Width in Feet 
(feet) : 100 :;: 300 : 500 : 700 : S500 : 500 : 500 
0 5092.5 5129.2 5144.6 5155.9 5166.5 5178.9 5193.9 
5 5038.6 5068.8 5083.9 5093.7 5101.6 5110.6 5119,2 
10 4982.9 5008.6 5023.6 5032.1 5036.8 5042.6 5047.8 
16 4921.1 4944.7 4958.0 4964.6 4967.2 4970.0 4972.8 
22 4859.4 4881.7 4892.3 4896.8 4898.6 4900.3 4901.8 
30 4777.7 4798.0 4805.3 4808.4 4810.1 4811.8 4812.2 
40 4675.9 4693.5 4698.9 4701.5 4702.9 4704.7 4704.5 
45 4625.2 4641.6 4646.4 4648.8 4650.1 4651.9 4651.6 
50 4572.6 4580.9 4583.0 4584.1 4594.4 4598.1 4598.6 
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Table 3 
DAMBRK ~ Prismatic 


oe, ” 
* t 


Pan wee 


Maximum Flow Depths Ne! 

~ y 

: a Maximum Depth in Feet a 

ee 

Distance : Manning Roughness Coefficient (n) fF 

From: 0.04 +: 0.04 : 0,06 : 0.04 : 0.06 : 0,08 : 0.12 Md 
Dam : Final Breach Width in Feet 


(miles) : 100 : 300 : 500 : 700 : 500 : 500 s 500 


0 62.5 99.2 114.6 125.9 136.5 148.9 163.9 
5 58.6 88.8 103.9 113.7 121.6 130.6 139.2 
10 52.9 78.6 93,6 102.1 106.8 112.6 117.8 
16 51.0 ee, 88.0 94.6 97.2 100.0 102.8 ue 
<i 
30 47.7 68.0 75.3 78.4 86.1 81.8 82.2 
re 
40 45.9 63.5 68.9 71.5 72.9 Wh.T 74.5 f: 
~~ 
4S 45,2 61.6 66.4 68.8 70.1 71.9 71.6 te 
50 42.6 50.9 53.0 54.1 64.4 68.1 68.6 
+ 
Ny 
: 
n 
fa 
be 
re. 
ry, 
‘, 
ll 
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Table 4 
DAMBRK — Prismatic 
Time to Peak Stage 


Time to Maximum Elevation in Hours 


Distance : Manning Roughness Coefficient (n) 
From : 0.04 : 0.04 : 0.04 =: 0.04 : 0.06 : 0.08 : 0.12 _ 
Dam : Final Breach Width in Feet 
(miles) : 100 : 300 : 500 : 700 : 500 : 500 : 500 
0 1.15 1.00 1.00 1.00 1.00 1.00 0.90 
5 1.45 1.15 1445 1.10 1.15 1.20 1.30 
10 1.70 1.35 1.30 1.20 1.40 1.50 1o75 
16 2.10 1.70 1,50 1.40 1.70 1.85 2.30 
22 2.45 2,00 ae oe 1.60 2.05 2.35 3.00 
30 2.95 2.35 2.10 1.95 2.55 3.10 4,00 
40 3.71 2.91 2.61 2.51 3.32 4.03 5.37 
50 4.52 3.35 2.99 2.83 4,09 5.07 6.89 
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FLOW SIM 1 RESULTS 


Computational instability was a major problem in applying FLOW SIM 1 to 
the prismatic case study. A solution could not be obtained with the original 
set of parameter values which included a 500-feet breach width, 1.0 hour 
breach time, and constant 0.04 Manning n. Numerous runs were made in an 
effort to obtain a solution with input data as close to these parameter values 
as possible. Alternative runs included various combinations of values for the 
following input. parameters: spatial step size (DELX), maximum hydraulic depth 
(HMAX) which is used by the program to determine the time step, Manning 
roughness coefficents (n), breach time, and breach width. A few runs were 
successful. Most runs terminated without obtaining a solution. The output 
from the unsuccessful runs usually included a message that the program 
terminated because of an instability in the calculations. 

Spatial step sizes of 1.0 mile, 0.5 mile, 0.2 mile, and 0.13 mile 
resulted in the same difficulties with computational instability. HMAX was 
varied from 40 feet to 400 feet with most runs using a value of 80 or 100 
feet. The effects of DELX and HMAX on stability were never clearly 
established. A DELX of 5,280 feet and HMAX of 80 feet seemed to perform as 
well as, and in some cases better than, other values of these parameters, 

Table 5 illustrates how obtaining a solution depends upon breach 
characteristics and n values, Several runs with alternative combinations of 
breach width, breach time, and constant Manning n are tabulated along with an 
indication of whether a solution was obtained or the run terminated due to 
instability. Input data other than the three parameters indicated were held 
constant for all the runs in the table. Values for DELX and HMAX of 5,280 
feet and 80 feet were used. 


A run with the same breach characteristics used in the Teton case study 
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(one-hour breach time, bottom width of 50 feet, and 0.66 sideslopes) and 
Manning n values of 0.040 also terminated without reaching a solution, 

The prismatic channel is representative of a steep-walled canyon. The 
cross-section side~slopes change abruptly from the canyon floor to the walls. 
An alternative run was made with the topwidths at the top of the canyon three 
times wider than the cross sections used in all the runs discussed above, 
This meant a much less drastic change in sideslope between the canyon floor 
and walls, With a breach width of 500 feet, breach time of 1 hour, and n- 
values of 0.04, DELX of 5,280 feet, and HMAX of 80 feet, the program still 
terminated due to an instability in the calculations. However, the stability 
characteristics of this widened prismatic channel were not investigated 
further. 

Table 5 summarizes the FLOW SIM 1 results for a breach width of 100 feet, 
breach time of 1.0 nour, and Manning n values of 0.06, The computer printout 
for this run is provided in Appendix B. The results of a DAMBRK run with 
these same parameters are presented in Table 7 for purposes of comparison. 
Figures 5 through 7 also provide a comparison of the FLOW SIM 1 and DAMBRK 
results for a breach width of 100 feet, breach time of 1.0 hour, and n values 
of 0.06. 

Appendix C provides a computer printout for one of the unsuccessful FLOW 
SIM i runs. This run used a DELX of 1056 feet, HMAX of 80 feet, breach width 
of 300 feet, breach time of 1.0 hour, and Manning n values of 0.060. The 
output includes a message that execution of the program was terminated because 


of an instability in the calculations, The output from the other unsuccessful 


runs is similar to the Appendix C printout. 
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Table 5 
FLOW SIM 1 -— Prismatic 
Computational Stability 


Breach Time : Breach Width : Manning n : Solution 
(hours) (feet) 
1 500 0.04 no 
{ 500 0.08 no 
1 500 0.12 no 
1 300 0.04 no 
1 300 0.06 no 
1 306 0.08 yes ¥ 
1 100 0.04 no 
1 100 0.06 yes 
1 100 0.08 yes 
2 500 0.04 no 
2 500 0.06 no 
2 500 0.08 yes 
2 300 0.04 no 
Z 300 0.06 no 
2 300 0.08 yes 
2 100 0.04 no 
2 100 0.06 yes 
2 100 0.08 yes 
5 500 0.04 no 
5 500 0.08 yes 
5 100 0.04 yes 
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Table 6 
FLOW SIM 1 - Prismatic Channel 
Results from a Successful Run 


Distance : : Maximum : Maximum : Time to 
Below : Peak : W.S. : Flow : Maximum 
Dam : Discharge : Elevation : Depth : Stage 


(miles) : (1000 cfs) : (feet msl) : (feet) +: (hours) 


1 995 5091.8 71.8 1322 
5 898 5047.2 67.2 1.53 
10 846 4989.9 59.9 2.06 
15 792 4938.0 58.0 2809 
20 748 4886.4 56.4 3.00 
25 711 4834.9 54.9 3.53 
30 677 4783.6 53.6 4.06 
35 648 4732.3 52.3 4.56 
40 621 4681.2 51.2 5.08 
45 295 4630.1 50.1 5.61 
50 575 4579.7 49,7 5.94 


Note: Breach width of 100 feet, breach time of 1.0 hour, 
Manning n values of 0.06, DELX of 5,280 feet, and HMAX 


of 80 feet. 
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Table 7 
DAMBRK - Prismatic Channel 
For Comparison with FLOW SIM 1 


Distance : : Maximum ;: Maximum : Time to 
Below : Peak ;: W.S. : Flow : Maximum 
Dam : Discharge : Elevatior : Depth > Stage 
(miles) : (cfs) : (feet msl) : (feet) : (hours) 
1 1,166 5095.5 75.5 1.30 
eS) 1,066 5050.5 70.5 1.55 
10 989 4992.6 62.6 1.90 
15 906 4940.0 60.0 2.38 
20 844 4888 .0 58.0 2.83 
25 797 4836.4 56.4 3.30 
30 760 4785.1 55.1 3.75 
35 730 4734.0 54.0 4.20 
40 705 4683.1 53.1 4.67 
45 682 4632.2 S252 5.13 
50 661 4580.7 50.7 5.69 


Note: Breach width of 100 feet, breach time of 1.0 hour, and 


Manning n values of 0.06. 
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FLOW SIM 2 RESULTS 


FLOW SIM 2 was also run with numerous combinations of values for the 
input parameters. Almost all runs terminated without obtaining a solution. 
The output from the unsuccessful runs usually included a message that the 
program terminated because of an instability in the calculations. Fssentially 
the same ranges of values for DELX, HMAX, breach time, breach width, and 
Manning roughness coefficient investigated wich FLOW SIM 1 were repeated with 
FLOW SIM 2. However, most of the combinations of input parameter values which 
ran successfully in FLOW SIM 1 terminated without solutions in FLOW SIM 2. 
The only successful runs of FLOW SIM 2 had a breach time of 5 hours. 

The printout for an unsuccessful FLOW SI 2 run is provided in Appendix 
D. The input data for this run includes a DELX of 1056 feet, HMAX of 80 feet, 
breach width of 100 feet, breach time of two hours, and Manning n values of 
0.06, The output from the other unsuccessful runs is similar to the Appendix 


D printout. 
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HEC-1 RESULTS 


Tables 8 and 9 and Figures 8 through 10 summarize the results of several 
HEC-1 runs with alternative sets of values for the parameters NSTPS and NMIN. 
The computational interval (NMIN) is not a critical parameter as long as it is 
small enough. Selection of a NMIN value is essentially a tradeoff between 
accuracy and computation costs. Since HEC-1 computer costs are modest, 
setting NMIN is not a problem. The number of steps (NSTPS) to be used in the 
storage routing is a somewhat arbitrarily determined parameter which can 
Significantly affect the computed hydrograph. The HEC-1 users manuai states 


that NSTPS is usually approximated by the following relationship. 
NSTPS = (reach length/average velocity)/time interval 


Computer runs are labeled by number in the following tables and figures, 
A value for NMIN of 4 minutes was used in runs 1 through 4 and 10 minutes in 
runs 5 through 7, NSTPS values used for each of 5 reaches were alternatively 
assumed to be aconstant 1, 5, and 10 for each reach. Runs were also nade in 
which the NSTPS values were estimated from the above relationship. For NMIN 
of 4 minutes, this resulted in NSTPS values of 4, 5, 7, 8, and 13 for the five 
reaches (run 4), For NMIN of 10 minutes, NSTPS values were estimated to be 2, 
2, 3, 3, and 5 (run 7). The five reaches are: miles 0 to 7.5, miles 7.5 to 
15, miles 15 to 25, miles 25 to 35, and miles 35 to 50, 

The results for NMIN values of 4 minutes and 10 minutes are very close. 
However, changing the NSTPS parameters significantly changes the resulting 
discharges and stages, Run 4 (NMIN = 4 minute and NSTPS values of 4, 5, 7, 8, 
13 for the five reaches) is probably the most realistic run, The computer 


printout for run 4 is provided in Appendix IE. 
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Table 8 a F 

HEC-1 - Prismatic Ni 

Peak Discharges and Water Surface Elevations ve 
. 
if 


Run ‘ 1 . 2 : 3 : 4 : 5 - 6 : 7 
NMIN : 4 : 4 : 4 $ 4 : 10 : 10 : 10 


NSTPS : 1 7 5 : 10 : varies : 1 5 : varies e 


. 
on Oe 
Distance : Maximum Discharges in 1000 cfs Ke 
(mile) ist 
0 3,911 3,911 3,911 3,911 3,911 3,911 3,911 % 
7.5 2,803 3,447 3,584 3,382 2,850 3,366 3,091 ny 
. 
15 2,159 3,162 3,428 3,110 2,182 3,338 2,641 ~ 
a 
25 1,542 2,596 3,107 2,704 1,553 2,798 2,180 ‘ 
35 1,183 2,140 2,617 2,344 1,193 2,317 1,811 nw 
6, ne 
€ 
50 799 1,507 1,823 1,856 803 1,576 1,428 : 
M3 
Miles :; Maximum Water Surface Elevations : 
‘ 
0 5301.7 5301.7 5301.7 5301.7 5301.7 5301.7 5301.7 — 
7.5 5094.7 5107.0 5109.5 5105.8 5095.6 5105.5 5100.4 . 
" 
15 4996.3 5014.7 5019.2 5013.8 4996.7 5017.7 5005.5 q 
a 
25 4895.5 4917.2 4926.3 4919.2 4895.8 4920.9 4909.2 * 
35 4786.7 4808.4 4817.6 4812.4 4787.0 4811.9 4801.5 a 
50 4650.8 4669.7 4677.1 4677.5 4651.0 4671.2 4667.9 a 
ia 4. 
‘ 
: 7 
§, 
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or 2 


Run 
NMIN 
NSTPS 


Distance 
(miles) 


Maximum Depth and Time to Maximum Depth 


a 


Table 9 


HEC-1 Prismatic 


in & rm 


Time to Maximum Depth in Hours 


0.02 


3 : 4 : 5 
4 : 4 : 10 
10 : varies : 1 


Maximum Depth in Feet 


154.5 150.8 139.6 
139.2 133.8 116.7 
146.3 139.2 115.8 
137.6 132.4 107.0 


147.1 147.5 121.0 


0.02 0,02 0.02 


1,20 1.20 1.17 


1.47 1.50 1.50 
1.87 1,93 2.00 
2.33 2.40 2.59 
3.13 3.13 3.33 
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0.02 
1.17 
1.50 
2.00 
2.90 


3.33 
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SMPDBK RESULTS 


The printout from the SMPDBK microcomputer program is presented as Table 
10. The results are plotted in Figures li through 13. Appendix F presents 


the results obtained from the manual computational procedure, 
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TABLE 10 


PROFILE OF CRESTS AND TIMES FOR 
SIMPLIFIED TETON DAM 
ON ALMOST PRISMATIC RIVER 


LONG OPTION 


DISTANCE PEAK PEAK TIME FLOOD TIME TIME FALL 
(MILES) FLOW DEPTH PEAK: STAGE FLOOD BELOW FS 
(CFS) (FT) (HRS) (FT) (HRS) (HRS) 
ie) 4015982 112.58 1 30 o1 1.5 
5 3234869 100.16 1.2 30 -S 2 
10 39211780 89.535 1.5 ete) «6 2.35 
25 2723214 81.77 2.2 39 1.3 3.2 
so 2216183 73.035 3.5 3 2.7 =) 


4x SUMMARY OF INFUT DATA #8 


NAME OF DAN 
NAME OF RIVER 


SIMPLIFIED TETON 
ALMOST PRISMATIC 


DAM CREST ELEVATION (FT.MSL) = 5301.7 
BREACH BOTTOM ELEV (FT. MSL) = S040 
FINAL BREACH WIDTH (FT) = S00 

VOLUME OF RESERVOIR (AC-FT) = 254192.872 
SURFACE AREA (ACRES) = 1940 

TIME OF FAILURE (MIN) = 60 
SPILLWAY/TURBINE FLOW (CFS) = 3580 


¥%a CROSS SECT*ON DATA 4% 


X-SEC NO. MILE NO. ELEV TOPWIDTH MANNING N) VALLEY HT FS 


1 G S030 oO 04 210 30 
5040 800 
5240 1400 
2 = 4980 =O 04 210 30 
4990 800 
5190 1400 
3 10 4930 fe) 04 210 30 Rs 
4940 1000 rg 
5140 1600 
4 26 4780 0 04 210 30 
4790 1000 
4990 1600 
5 SO 4530 0 04 210 30 
4540 1000 


4740 1600 
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HEC DIMENSIONLESS GRAPHS PROCEDURE RESULTS 


The HEC dimensionless graph procedure is based on the assumption of an 
instantaneous complete removal of the dam. Thus, the breach characteristics 
(500 feet width and 1.0 hour breach time) used in the other models could not 
be incorporated in the dimensionless graph procedure. The prismatic channel 
characteristics are the same as used in the other models. The required 


computations for the procedure are presented in Table 12 and the results are 


summarized in Table 11. 
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Table 11 
HEC Dimensionless Graphs - Prismatic 
Summary of Results 


Distance : Maximum : Time to 
From : Flow : Maximum 
Dam : Depth : Depth 

(miles) : (feet) 3 (hours) 
0 119.5 0.00 
5 101.6 0.43 
10 95.8 0.68 
15 91.3 0.79 
20 87.4 1,07 
25 84.5 1.18 
30 81.5 1.36 
40 IS 1.75 
50 72.8 Ze22 


Table 12 
Computations for HEC Dimensions Graphs Procedure 


So = 10ft/mile = 0.001894 
Yo = 5310.7 - 5030 = 271.7 feet 


Lo = Yo/So = 271.7 £t/0.001894 = 143,458 ft = 27,17 mile 


Head : Channel : Pool : Reservoir Volume 
Point : (feet) :; Elevation : Elevation : (acre-feet) : (£t3) 
] 271.7 5030 5301.7 251,300 1.095 x 1010 
2 120.8 5030 5150.8 50,000 2.178 X 109 


M= (log (Volo / Vol,) / log (¥/Y1)) - 2 
M = (log (50,000/251,300) / log (120.8/271.7)) - 2 
M = -0.008 = 0 


ze S “| 
C = (Vol;)(3)(M+1) | —~ 9. 
(Vol, )(3)(M+1) a7 
Y} 

C = (1.095 X 1019)(3)(1) sore on = 842.6 
oye ey, 

C = (2.178 X 109)(3)(1) | 2.001894 | = g48.1 
120.82 

use C = 845.4 


C= C(Y¥M-1) = 845.4 (271.771) = 3.11 
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Table 12 Continued 


Computations for HEC Dimensionless Graphs Procedure 


Bp = 1,585 ft 

Yo = 271.7 ft 

Ko = Bo * Yo = (1,585)(271.7) = 430,645 fe? 
Py = Bo + 2¥o = 1,585 + 2(271.7) = 2128 ft 
Ry = Ap/Po = 430645/2128 = 202.3 ft 


Vo = 1:486 Ro2/3 Sol/2 = 1.486 (202,3)2/3(.001894)1/2 
n n 


55.7 fps 


< 
° 
N 


Fe Vs 55.7 = .596 = .60 
(go Yo —-('32.2(271.7) 


Lot Yo/So = 271.7/.001894 = 143,453 ft = 27.2 mile 


2575 seconds 
0.7153 hours 


T, = To/Vo = 143,453/55.7 


Fo = e 60 


Lo = 27.17 mile 


To = 42.9 hours 
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Table 12 Continued My 
Computations for HEC Dimensionless Graphs Procedure a 


Me 
Maximum Flood Depth a 


Time of Occurrence of Maximum Flood Level ad 


X : X : Fo=0.5 : Fo=1.0 ; Fo=#0.6 : tw (hours 


0 0 0 0 0 0 

5 0.184 0.5 0.7 0.6 0.43 
10 0.368 0.8 Led 0.95 0.68 aay 
Yi 
15 0.552 1.0 12 1.1 0.79 , 
ny 
20 0.736 1.3 ligt 1.5 1.07 ph 
25 0.920 1,5 i.8 1.65 1.18 a 
a 
30 1.104 1.8 2.0 1.9 1.36 4 
S\9 
40. 1.472 2.3 2.6 2.45 1.75 a 
50 1.840 3.0 3.2 3.1 199 . 
ae 
a 
40 oi 
ie 


| SS 


- 
~- 


as toe 


SIttW UL WEG MOLDG JIUCISIC 


0S Se OF sf ae S2 02 St Ol S 0 


Of 


Ob 


0s 


09 


0d 


4} 


0g 


2994 Uj YyWdad WNuxeYN 


06 


00l 


Olt 


suzdaq MOL4 4eed 
DtyeusStuy - sydeug ssaluolsuaieg 


pl aunbe4 . oe 


Sa(iw ut weg Molag aoueqsey 


| | 5 
| a 2 
‘ } = 
| | = 
. 
H 3 
| | m 
2 nN 
: “ee Psy +t 
| | S 
s 
| | a 
| he | 3 
i made ujdag wnwexey 02 auey | ~ | 
DL 7ZeWStug - Ssydeug ssajuoLsuautg | | 
| Gt aanbey | 


TR66 RESULTS 


The SCS Technical Release No. 66 simplified dam breach routing procedure 
results are presented in Table 13 and Figures 16 and 17. The calculations are 
Summarized in Table 14. This manual computation procedure provides peak 
discharges and depths at selected locations. The peak discharge at the dam is 
determined from a relationship between peak discharge and reservoir water 
depth which was developed based on data from actual past dam failures, The 
peak discharge of 1,931,000 cfs computed with the TR66 procedure is about half 
the magnitude of the peak outflow from the 500-feet wide breach computed by 
the other models. Consequently, an additional simulation was made in which 
the 3,841,000 cfs peak discharge at the dam computed with the DAMBRK model was 
routed downstream with the TR66 valley routing procedure. The results using 
TR66 for both breach outflow and valley routing are labeled TR66 in the tables 
and figures. The label TR66 (DAMBRK) is used to signify that the DAMBRK peak 


breach outflow is combined with TR66 valley routing. 
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Table 13 
TR66 - Prismatic 
Peak Discharges and Flow Depths 


Distance : TR66 : TR66_(DAMBRK) 
From Dam : Discharge : Depth : Discharge : Depth 
(miles) : (cfs) : (feet) : (cfs) : (feet) 
0) 1,931 84.5 3,841 123.0 
5 1,506 73.5 2,996 107.9 
19 1,313 62.3 2,612 89.0 
15 1,130 57.5 2,247 82.0 
25 834 48.7 1,659 69.0 
35 715 44,2 1,421 64.2 
50 525 37.5 1,045 55.0 
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tf 
Table 14a 0. 
TR66 Calculation Summary 
Depth Versus Area and Volume Relationships 

Depth : Miles Downstream From Dam q 
(feet) : Q : 5 : 10 : 25 : 50 = 

Cross-Sectional Area (ft2) M3 

0 0 0 0 0 0 t 
10 4,900 4,000 5,060 5,000 5,000 a 
60 47,750 47,750 58,750 58,750 58,750 : ; 
110 99, G00 99 ,000 120,000 120,000 120,000 : 
160 157,750 157,750 188,750 188,750 188,750 4 
Incremental Volume (million ft) .: 

0 0 0 0 0 i 
10 105.6 118.8 396.0 660 a 

60 1,261 1,406 4,653 7,155 ; 
110 2,614 2,891 9,504 15,840 i 

160 4,165 4,574 14,949 24,915 ‘ 
Total Volume (million ft3) 2 
0 0 0 0 0 0 a 
10 0 105.6 224.4 620.4 1,280.4 7 
60 0 1,261 2,667 7,320 15,075 “" 

110 0 2,614 5,505 15,009 30,849 : 
160 0 4,165 8,739 23,688 48,603 a 
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Table 14b 
TR66 Calculation Summary 
Depth Versus Discharge Relationships 


Q = 1.49 ar2/3 sl/2 = As40 AR2/3 (,0018939)1/2 = 1.6211 AR2/3 
n 


° 


mile 0.0 to 5.0 


: Wetted : Hydraulic : 
Depth : Topwidth : Area (A) : Perimeter (P) : Radius (R) : Discharge (Q) 
(ft) (ft) (£t2) (ft) (ft) (cfs) 
0 0 6) 0 6) 0 
10 800 4000 800.25 4,998 18,955 
60 950 47,750 980.53 48.698 1,032,268 
110 1100 99,000 1160.58 85.302 3,109,943 
160 1250 157,750 1340.85 117.65 6,140,045 
mile 10.0 ta 50.0 
: Wetted : Hydraulic : 
Depth : Topwidth : Area (A) : Perimeter (P) : Radius (R) : Discharge (Q) 
(ft) (ft) (£2) (ft) (ft) (cfs) 
0 0 6) ) 0 6) 
10 1000 5000 1000.2 4.999 23,697 
60 1150 58,750 1180.5 49.767 1,288,584 
i10 1300 120,000 1360.8 88.183 3,854,036 
16C 1450 388,750 1541.0 122.485 7,546,630 


Table 14¢ 


TR66 Calcvlation Summary 
Outflow Hydrograph 
Qnax = 65 Hl-85 = 65 (261,7)1-85 = 1,931,300 cfs 
3 
Q.2 = BA = (32.2) (74,000) 
i T 1,023.5 


Q = 3,570,520 


Since flow is subcritical, the curvilinear hydrograph is used. 
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Table 14d 
TR66 Calculation Summary 
Maximum Discharge and Depth 


n N 
(N-n)S_ log Qi -n D> log Q 
i=1 ien+t+l 
n- N 
(N-n) >_ log Sjeicn >_ log Sj,i 
i=l i=an+l 


n th 
_ log Q - m>_log Sjri 
K = log-l | isl isl 


n 


a 

— — = 1,931,300 cfs 
2 Kva K (10,947 million ft3)@ 
* 


* x 
Ky (1 - Q,)65/3 - @) from figure 2 of TR66 


vw 
ul 


x 
Qo and Q. from ES-212 of TR66 


d from curve of Q versus d 
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Table l4e P 

TR66 Calculation Summary ; 

Maximum Discharge and Depth F 

Na 

Vv 

* * * * ve 

Location: m : K : Ky : Q : KF : Q : Q : od ‘ 

(miles) (cfs) (feet) i" 

i 

0 ES a 5 os = - 1,931,300 84.5 y 
5 1.563 13.84 0.06783 0.770 0.0582 0.780 1,506,000 73.5 
10 1.564 5.306 0.1752 0.665 0.157 0.680 1,313,000 62.3 
15 1.569 2.492 0.3562 0.540 0.330 0.585 1,130,000 57.5 
25 1.573 1.015 0.8425 0.450 0.946 0.432 834,000 48.7 
35 1.593 .4956 1.433 0.360 1.387 0.370 715,000 44.2 
50 1.576 .3175 2.619 0.265 2.547 0.272 525,000 37.5 
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COMPARISON OF MODEL RESULTS 


A comparison of the results obtained using DAMBRK, HEC-1, SMPDBK, TR66, 
and the HEC dimensionless graphs procedure are presented in Tables 15 through 
18 and Figures 18 through 20. Both HEC-1 runs 1 and 4 are shown. HEC-1 run 1 
used a NSTPS value of 1 for all reaches, whereas run 4 used NSTPS values of 4, 
5, 7, 8, and 13 for the five reaches. The HEC dimensionless graph procedure 
assumes an instantaneous complete rem.val of the dam. The other models 
reflect a breach width of 500 feet and breach time of 1.0 hour. Solutions 
could not be obtained with these breach parameters using FLOW SIM 1 and FLOW 
SIM 2, Runs with these two models terminated with a message that an 


instability in the calculations had occurred. 
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Table 15 
Comparison of Models - Prismatic 
Peak Discharges 


Distance : Maximum Discharge in 1000 cfs 


Fron Dam : DAMBRK : 4HEC-1 :  HEC-1 =: SMPDBK : TR66 :;: TR66 

(miles) : : (Run-1) +: (Run-4) : 7 : (DAMBRK) 
0 3,841 3,911 3,911 4,016 1,931 3,841 
5 3,468 3,172 3,508 35235 1,506 2,996 
10 3,220 2,588 3,291 3212 1,313 2,612 
15 2,925 2,159 3,110 3,049 1,130 2,247 
20 2,529 1,815 2,910 2,880 985 1,950 
25 2,380 1,542 2,704 2,723 834 1,659 
30 2,135 1,340 2,520 2,620 780 1,540 
35 1,936 1,183 2,344 2,520 715 1,421 
40 1,777 1,040 2,180 2,410 660 1,300 
45 1,770 920 2,020 2,310 600 1,175 
50 1,567 799 1,856 2,216 §25 1,044 

: 

i: 
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Table 16 
Comparison of Models ~ Prismatic 
Peak Water Surface Elevations 


Distance : Maximum Water Surface Elevation in Feet ms]. 


From Dam : DAMBRK : HEC-1 : HEC-1 : SMPDBK : Dimensionless : TR66 : TR66 


30 4583.0 4650.8 4677.5 4603.0 4602.8 4567.5 4584 4 


(miles) ; : (Run-1) : (Run-4) ; : Graphs : : CDAMBRK) 

9 5144.6 - - 5142.6 5149.5 5114.5 5153.0 

5 5083.9 - - 5080.2 5081.6 5053.5 5087.9 

10 5023.6 5061.9 5075.1 5019.5 5025.8 4992.3 5019.0 

15 4968.9 4996.3 5013.8 4966.9 4971.3 4937.5 4962.0 

20 4913.9 4945.5 4965.8 4914.0 4917.7 4883.0 4905.3 

! 25 4859.6 4895.5 4919.2 4861.8 4864.5 4828.7 4849.0 

: 30 4805.3 4840.8 4866.0 4809.9 4811.8 4776.4 4796.4 

35 4751.8 4786.7 4812.4 4758.3 4758.4 4724.2 4744 ,2 

40 4698.9 4741,1 4767.8 4706.7 4708.1 4672.1 4691.0 

45 4646.5 4696.0 4722.6 4654.5 4654.U 4620.0 4638.0 
\ 
\ 
| 


Table 17 
Comparison of Models - Prismatic 
Maximum Flow Depths 


ed 


Distance ; Maximum Flow Depth in Feet 


From Dam ; DAMBRK : HEC~l : HEC-1 : SMPDBK : Dimensionless : TR66 : TR66 

(miles) : : (Run 1) : (Run 4) ; ; Graphs : : CDAMBRR) 
0 114.6 - - 112.6 119.5 84.5 123.0 
5 103.9 - - 100.2 101.6 73.5 107.9 
10 93.6 131.9 145.1 89.5 95.3 62.3 89.0 
15 88.9 116.3 133.8 86.9 91.3 57.5 82.0 
20 83.9 115.5 135.8 84.0 87.7 53.0 75.3 
25 79.6 115.5 139.2 81.8 84.5 48.7 69.0 
30 Tae3 110.8 136.0 79.9 81.8 46.4 66.4 
35 71.8 106.7 132.4 78.3 78.4 44.2 64.2 
40 68.9 111.1 137.8 76.7 78.1 42.1 61.0 
45 66.5 116.0 142.6 74,5 74.0 40.9 58.0 
50 53.0 120.8 147.5 73.0 72.8 37.5 54,4 

- 

R 

: 
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Table 18 e 
Comparison of Models — Prismatic 3 
Time to Peak Stage t 
. 
Distance: —s—sé<‘__ Time to Maximum Depth in Hours 4 
T'rom Dam : DAMBRK : HEC-1 : HEC-1 : SMPDBK : Dimensionless 
(miles) : : (Run 1) : (Run 4) : : Graphs i 
4 
0 1.00 - ~ 1.0 0.00 : 
5 1.15 - - 1,2 0.43 = 
10 1.30 1,24 1.30 1.5 0.68 ‘ 
15 1.46 1.47 1,50 1.7 0.79 ; 
20 1.65 1.70 1.72 1.97 1.07 : 
25 1.85 1.93 1.93 2.2 1,18 S 
30 2.10 2.20 2.17 2.45 1.36 7 
35 2.35 2.47 2.40 2.7 1,56 
40 2.61 2.75 2.76 2.98 1.75 
45 2.88 3.04 3.04 3.24 1.98 
50 2.99 3.33 3.13 3.9 2.22 


56 


a : ‘ : ‘ 3 i lee : , 4 poe «. ] oe -; ‘ * «. , a . woe 

A iy ge <1 he le Oe em 5. gs a See ee? rs ae he id i a Py od wg 0 5 f Aft br a a ee: 

pe wraeges NN lag Nd | | ot witeef ate: fos ee 3 re, = ma x ; So, x9 a ee aS q e! Ae bees . x : “= | 
beth cee oe ener ios Pong wars I ES, be Teall eh 30 FP. Pa Sa 53 - Le eo eo me Y Ae a? oA GRU + ko 76". “all : on 


SOLiW Ut WEQ MOLag adueysSIG 


0S Sb Op Se o€ S2 C2 Sl ot ) 0 


vs ie BONA ah oat | Rene he NE, CE gman A te a eee PET Ate eh ne 
e as ea Sew I Ee Te » ek, Say Ye cea 
met oes Pee eld ‘enw Sconce iokanle Aap en is sosets pee aa Ralibaat ivan, wag ae Faas 
: Sao ede PELENaSeL be Cie an cat (a ard [Ca a aA Sain Pets Gar add eyes sand Sea 4 
fain tnty to wel hag eet Cee bye Sa Eee ae it te ee eel ere eee Weise 
Beare oye tan eee 7-0 is orem Reece tet eet ee ees ees 
it + ae ‘ ee e4 hed % boas go ate se rs cee ie rr ry 
7 ides oar Ce ee eee ae ero pe cts ee tee ag rer 
; aaah ie ahaa eae Gece. Se eal (iad, aoe bs RS eer haps 


W) £-9274 


H "“ : 
7 halen 
: M a 
Lz 
J —s Sy 
we 
aly: Es 
1 


00S‘ I 


¢ wed Le 790! WMIGIT z 
| Med XYGWYO) 2 > 
SS : 
Le ie | ; Z ~ 
2 Sel ee nee CPC eS 0002 a 
~~ MSALEN FITVIAKA & 
— APP. 3 
3 
00S'2 2 
Oo 
° 


o00°€ 


* ‘ 
sabueyrdsiq eed | 
i ; DLzeWSLug - UOSLUediIO) LapoW | | i 
Bl aundt4 | : 

| 
1 


nn et 600' 2 


SaliW UL Weg MOLag azuveysIG 


0 

syzdaq MC. 4 4Rag 01 
JLJeWSlug - UOSLUedWO) [apoy 

6L aunbey : G2 

Of 

bas |e Ov 

0S 

(WO LY. eet ePdd. |. 8? 09 


MOK MYPAUIY 99 HL 


ra ‘ | 
HdASYD SSYINOISNAW, oan 7 | a 
; ; ! 


IY 
SALSN FIGHISVA) 6-27 


3984 Ul Yidad whuLKey 


4 


ieaks! een] SAR RAGRAS: sogseee | inet etese | \eeeetes) Vanes Oe 


” al 


ae 


SO,tW UL weg MOL ag a2uUeqSIG 


0S St Op SC ag 2 a2 St Ol S 0 
) 
| 
| he ee Ds, woe ne 5‘0 
enc lts atetoe er ee ae eee o'k 
! | 
| , 
| 
SHIOND SSTINGISNINIT 2.31 ~.. ! Peep be S*t 
f 
| | | 
| | | . 
-CSALSN F7AMIUA) [-~22H' ae 


SUPOW UL YOLJEAD[ 5 Wu Key] OL Wty 


| 
| 
| 


yjdaq wnwixey 07 owt) 
Diyewstug - uostseduios) Lapow 
O2 aunbr4 


| 
| 
| 


be Nl A ar Sh ae 


Appendix A 
DAMBRK Base Run 
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Appendix F 


SMPDBK Manual Computations 


2 ae a 
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Manual Version of SMPDBK 
Hypothetical Prismatic Channel Case Study 
Summary of Calculations 


Reference: Jonathan N. Wetmore and Danny L. Fread, "The NWS Simplified 
Dam-Break Model Executive Brief." 


INPUT DATA 


> 
iT} 


1940 acres 


on 
" 


5301.7 - 5040 = feet 


& 
" 


500 feet 


te = 1.0 hour 


> 
o 
ul 


3,580 cfs 


Vol = 251,300 ae - ft = 1.09466 x 101Y ¢43 


mile 0.0 


elev : depth : topwidth 


5930 0 0 

5040 10 800 

5240 210 1400 
mile 10.0 


elev : depth ; topwidth 


4930 0 0 

4940 10 1°00 

5140 210 1600 
mile 50.0 


elev : depth : topwidth 


4530 0 0 
4540 10 1000 


4740 210 1600 


F3 


mile 5.0 


elev : depth : topwidth 


4980 0 0 

4990 10 800 

5190 210 1400 
mile 25.0 


elev _: depth : topwidth 


4780 0 0 
4790 10 1000 


4990 210 1600 


ee ES a aa 


CHANNEL DESCRIPTION 


(Zlog hi)(Zlog Bi) 


m = (Clog hi)(1cg Bi)) - I 
> (log hi)? - (Glog hi)? 
I 


log k = Lhog Bi _ n 2408 hi 


For Miles 0 to 5 


(3.322)(6.049) 
ip 209 ee RCT 


6.393 - _(3,322)¢ 
2 


log k = o.0ee - 0.18466 3.322 = 2.7177 


k = 522 


B= kh™ = $22 70-1847 


For Miles 10 to 50 


B = 675 h0+1577 


F4 


MAXIMUM BREACH OUTFLOW (Somax) 


Ca 23.4 As = 23.4 020) = 90.792 


By 
a ae 
Qbmax = % = 3.1 Br | -—~ 
le fH 
|___ 90.792 
= 3,590 + 3. 100) 90753 
261.7 
4 21869 0.05 ok 1.49110. | 
n (m1)>/3 0.04 | 5280 | 


b sm + 5/3 = 0.1847 + 5/3 = 1.851 


11/1.851 
h = os i | 4 4,016,000 : 
cael "638.2 2 
Check for Submergence 
m a ee . 90.792__| 
weir ltr + © | 1.0% 90.792 | 
iY | | 7261.7 | 


Submergence correction not required. 


i 522 Mie 
| (ager) 573 


= 112.7 ft 
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DOWNSTREAM ROUTING 


mile 5 
xX 26,400 
FC 0.605 
\* 1.326 
XC 73,809 

X/XC 0.358 


QP/Quax «0+ 922 
QP 3,703 


Mile 


29 
50 


ee en es 


10 
52,800 
0.605 
1.324 
72,700 
0.726 
0.851 
3,418 


4,016,000 
3,703,000 
3,418,000 
2,827,000 


2,157,000 


A ee ee a 


25 
132,000 
0.606 
1.319 
69,500 
1.899 
0.704 


29827 


: Peak Discharge (cfs) 


50 
264 ,000 

0.607 

1.310 


64,000 


- 


LSSESS. 


ares 


—_ 
—_ 


